Full quantum distribution of contrast in interference experiments between interacting 

one dimensional Bose liquids 
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We analyze interference experiments for a pair of independent one dimensional condensates of 
interacting bosonic atoms at zero temperature. We show that the distribution function of fringe am- 
plitudes contains non-trivial information about non-local correlations within individual condensates 
and can be calculated explicitly using methods of conformal field theory. We point out interesting 
relations between these distribution functions, the partition function for a quantum impurity in a 
one-dimensional Luttinger liquid, and transfer matrices of conformal field theories. We demonstrate 
the connection between interference experiments in cold atoms and a variety of statistical models 
ranging from stochastic growth models to two dimensional quantum gravity. Such connection can 
be used to design a quantum simulator of unusual two-dimensional models described by nonunitary 
conformal field theories with negative central charges. 



I. INTRODUCTION 

The hallmark of a Bose-Einstein condensate (BEC) is 
the existence of a well defined macroscopic phase. In- 
deed experiments with large condensates display robust 
matter wave interference with negligible fluctuations in 
the fringe contrast 1 . From this point of view, large three 
dimensional condensates may be thought of as classical 
objects. However, there is a continuous range of possibili- 
ties intermediate between perfect condensates on one side 
and systems that do not display an interference pattern, 
like high temperature thermal gases, on the other—. For 
example, one-dimensional interacting bosons display in- 
terference patterns with a reduced contrast and with non 
negligible shot to shot fluctuations of the fringe contrast. 
In this paper we analyze the full distribution associated 
with this quantum noise, and discuss what it can tell us 
about the underlying strongly correlated state. 

By virtue of its direct connection to the concepts of 
quantum measurement, study of quantum noise has deep- 
ened our understanding in a variety of areas. Under- 
standing the noise in photo detection 4 prompted the cre- 
ation of nonclassical states of light and led to the de- 
velopment of quantum optics. In mesoscopic electron 
systems, current fluctuations contain information that is 
not available in simple transport measurements. For ex- 
ample, they can be used to distinguish electrical resis- 
tance due to diffusive scattering from the resistance due 
to point contact tunneling (see Ref. [j| for a review) . Sin- 
gle atom detectors have been used recently to perform 
Hanburry-Brown-Twiss experiments with cold atoms§i 
Finally, analysis of noise correlations in time of flight ex- 
periments of ultra cold atoms has been proposed— and 
tested*^, promising a powerful new technique to access 
many-body correlations in such systems. In these appli- 
cations, the quantities of direct interest are contained in 
the first few moments of the noise distribution, such as 
the noise power spectrum. However deeper insights into 
quantum systems may be gained by obtaining the full 



statistics of the fluctuations. 

One of the central problems in the field of ultra cold 
atoms is finding new ways of characterizing many-body 
quantum states. In this paper we demonstrate that anal- 
ysis of the distribution function of contrast in interfer- 
ence experiments between interacting one dimensional 
Bose liquids provides a novel probe of non-local corre- 
lations and entanglement present in these systems. It is 
known that the average amplitude of interference fringes 
can be used to extract two point correlation functions in 
fluctuating condensates 2 . This idea has been successfully 
employed by Hadzibabic et al. to measure the Kosterlitz- 
Thouless transition in two dimensional systems^. Inter- 
ference experiments, however, contain more information 
than the average value of the contrast. Each observation 
of the interference pattern is a classical measurement of 
a quantum mechanical state, so the result of each indi- 
vidual measurement will be different from the average 
value. As we discuss below, higher moments of the dis- 
tribution function of interference amplitudes correspond 
to high order correlation functions. Hence the knowledge 
of the entire distribution function reveals global proper- 
ties of the system that depend on non-local correlation 
functions of arbitrarily high order. So far, the discus- 
sion of full counting statistics has been limited mainly to 
systems of non interacting particles 1 —. The main result 
of this paper, by contrast, is an expression for the full 
distribution of amplitudes of interference fringes arising 
from one dimensional interacting Bose liquids, that can 
be described by a Luttinger liquid. 

The approach used in this paper for analyzing interfer- 
ence experiments can be generalized to a variety of other 
measurements in cold atom systems. For example, one 
can analyze fluctuations in particle number in systems 
with pairingii and fluctuations in magnetization in Mott 
states of atoms with magnetic exchange interactions. In 
both cases distribution functions will contain non-trivial 
information about underlying many-body states. One 
can also consider time dependent phenomena such as evo- 
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lution of phase coherence between a pair of condensates 
coupled by a finite tunneling amplitude 3 . In this pa- 
per we focus on the distribution function of contrast in 
interference experiments between independent Bose liq- 
uids in their ground states. Already in this case we find 
a non-trivial evolution of the distribution function from 
the non-fluctuating perfect contrast for the case of non- 
interacting atoms to the Poissonian distribution of con- 
trast for atoms in the regime of infinitely strong repul- 
sion. We think that our work constitutes one of the first 
steps in the direction of developing a new method for 
characterizing interacting many-body quantum states of 
cold atoms using full distribution function of some exten- 
sive quantum operator, such as the total magnetization 
or the interference amplitude operator defined in equa- 
tion JIJ. Different kinds of cold atom systems and ex- 
perimental probes can be studied following this general 
approach. 

We start our analysis by establishing the relation be- 
tween the probability distribution function of the inter- 
ference amplitude and the partition function of a bound- 
ary sine-Gordon model. Using methods of conformal 
field theory, we reduce this problem to that of finding a 
spectral determinant of a simple one-dimensional, single- 
particle Schrodinger equation. We solve this problem nu- 
merically, as well as analytically using the WKB approx- 
imation, to obtain the desired fringe distribution for any 
value of the Luttinger parameter. 




FIG. 1: Schematic view of the possible experimental setup, 
which produces interference pattern between two independent 
ID condensates. 

It is interesting to point out that one can take the 
alternative point of view on the relation between inter- 
ference experiments with BEC and the quantum impu- 
rity problem. By measuring the distribution function of 
the interference amplitudes experimentally one obtains 
the full partition function of the boundary Sine-Gordon 
problem (see eq. © below). This model describes a 
range of interesting problems such as an impurity in a 



Luttinger liquid^*, the asymmetric Kondo modeU^, and 
the tunneling of a particle in the presence of dissipation 
within a Caldeira-Leggett approach^. As we discuss be- 
low, interference experiments can be used to obtain the 
partition function of these systems not only in the ground 
state but also in the non-equilibrium regimes (e.g. in the 
presence of a finite voltage). Hence interference exper- 
iments can be considered as a quantum solver of these 
non-trivial many-body problems. 

Another surprising implication of our analysis is 
that interference experiments with one dimensional cold 
atoms can be used as a quantum simulator of several 
fundamental problems in physics. This connection relies 
on the fact that many models of systems with critical 
behavior in two-dimensional statistical mechanics, one 
dimensional field theories and many-body quantum sys- 
tems can be described by continuum theories with con- 
formal invariance. Basic ingredients of any conformal 
field theory are the central charge and conformal dimen- 
sions of operators. One dimensional quantum systems 
of interacting particles have positive central charges and 
correspond to the so-called unitary class. On the other 
hand, there is a non-unitary class of models which have 
negative central charges. Such models appear in con- 
texts as different as two dimensional quantum gravity 
and stochastic growth models and they have very un- 
usual properties (see e.g. Ref. 17]). In this paper we 
demonstrate that interference experiments with one di- 
mensional condensates can be used to analyze models 
in nonunitary universality class by virtue of the exact 
mathematical relation between the distribution function 
of interference amplitude and the so-called Q-operators 
of conformal field theories. In Fig.0]below we show a dia- 
gram which illustrates several models that can be studied 
using interference experiments. 



II. FULL QUANTUM DISTRIBUTION OF THE 
INTERFERENCE CONTRAST 

The setup for the interference experiments we consider 
is sketched in Fig.^ Two independent quasi-condensates 
are allowed to expand in the transverse direction. After 
sufficient time of expansion, the integrated density pro- 
file is measured by the imaging beam which is sent at 
an angle 9 to the condensate axis. This setup is quite 
common for the interference experiments in cold atoms 
and was already realized by several groups pdHT^ |. 

Everywhere in this paper we consider the two con- 
densates to be identical, although our analysis can be 
generalized to the case of different condensates (see also 
Ref. We assume that before the expansion, atoms 

are confined to the lowest transverse channels of their 
respective traps and that the optical imaging length L 
(which is smaller than or equal to the size of the sys- 
tem in the axial direction) is much larger than the co- 
herence length of the condensates. This allows us to use 
an effective Luttinger liquid description of the interacting 
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bosonsSi. The operator corresponding to the interference 
signal of the two condensates^ is given by 



dza[(z) 02(2) e Vj 



(1) 



Here a\ and 02 are the bosonic operators in the two sys- 
tems before the expansion, and the integrals are taken 
along the condensates. The seeming winding of the rel- 
ative phase between the two systems, described by the 
exponential term in equation Q), can either come from 
the measurement process itself or from the actual motion 
of the condensates 2 . If the condensates are at rest, we 
have p = (md/ hi) tan 9, where m is the atom's mass, d 
is the separation between the two condensates and t is 
the time when the measurement was done after the free 
expansion started. With some abuse of terminology we 
will call p the relative momentum. When condensates 
one and two are independent, one finds that the expec- 
tation value of (Ap) vanishes, however (|^4p| 2 ) is finite. 
This means that individual measurements show a finite 
amplitude of interference fringes, however their phase is 
completely unpredictable. Higher moments of the inter- 
ference fringe amplitude are given by 2 - 
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Al n Z<fJ, where.! 



o = y/cp^L^, (2) 



where C is a constant of order unity, p is the particle 
density in each condensate, £/j is the short range cut- 
off equal to the healing length, and K is the Luttinger 
parameter describing the interaction strength. In this 
letter we assume that K > 1, which is always the case 
for bosons with S- function type repulsive interactions 38 . 
Coefficients Zm in Eq. (J2J) are given by^i: 



n i<J 2sin(^)n fc<; 2sin(^) 



i2p ^2 ■ (ui —Vi) 



IL, fc 2sin (- 



l/K 



(3) 



where p is the relative momentum measured in units of 
2tt/L: p = pL/2n. 

Coefficients originally appeared in the grand 

canonical partition function of a neutral two-component 
Coulomb gas on a circle 



Z p (K, x) 



^ (n!) 2 



(4) 



Here x is the fugacity of Coulomb charges and Zm de- 
scribe contributions from configurations with 2n charges 
(i.e. canonical partition functions). The partition func- 
tion 0} with K > 1 and p being half-integer describes 
several problems in statistical physics (see Ref. j2^ and 
references therein). In particular, it describes an impu- 
rity in a one-dimensional interacting electron liquid. At 



low energies this problem is described by a Luttinger liq- 
uid (LL) with an additional local non-linear term due to 
backscattering from the impurity: 



S = — / dy dr [(d T 
1 J -co Jo 



C») 2 + (dy<t>f 



2g / dr cos 
Jo 



2tt0(O,t) + |t 



,(5) 



where g is the amplitude of backscattering on the im- 
purity and 4>(x) is a bosonic phase field associated with 
the electron field operator ip(x) ~ e 1 ^ . In the bosonized 
form the electron-electron interaction becomes quadratic 
and is effectively described by the Luttinger parameter 
K. Perturbative expansion of the corresponding parti- 
tion function in powers of g produces the series (0} with 
the fugacity given by x = <7/3(27t//3k) 1 / 2K . Here k is a 
non-universal renormalization factor, which sets the scale 
for the long distance asymptotics of the correlation func- 
tions: (exp[2£(0(i,r) - 0(0,0))]) - (nVx 2 + t 2 )-^ k . 
Finally, the single impurity Kondo model is related to 
Z p (K,x) as welU£. 

It is easy to understand the origin of the relation be- 
tween interference experiments and a quantum impurity 
problem. Moments of fringe amplitudes are determined 
by high order correlation functions computed at the same 
time but in different points in space. On the other hand, 
expansion of the partition function for a quantum impu- 
rity contains correlation functions computed at the same 
spatial point but at different times. Lorentz invariance of 
the LL ensures that the two are the same. Note that the 
analogue of the finite imaging angle 9 in the interference 
experiments is a finite voltage in the quantum impurity 
problem 2 ^. This analogy can be also understood from the 
interchanged roles of space and time in the two systems. 

When describing interference experiments it is con- 
venient to define the normalized amplitude of interfer- 
ence fringes a — A^/Aq. From eq. J2} we find that 

(a n ) = Z^ , so by measuring the distribution function 
Wp(a) experimentally, we get direct access to the parti- 
tion function (gj. We point out that W p (a) can be used 
to compute all moments of Ap, and therefore contains 
information about high order correlation functions of the 
interacting Bose liquids. 

Using Taylor expansion of the modified Bessel function 
as well as Eq. (0} and the fact that (a n ) — Z% n we find 

Z p (K,x)= W p (a) I (2xVa) da. (6) 
Jo 

Inverting Eq. 10 we can express the probability W p (a) 
through the partition function Z p (K,x). Noting that 
Io(ix) — Jo(x) and using the completeness relation for 
Bessel functions, J Q Jo(Ax) Jo(Xy)\x\\d\ = S(\x\ — \y\), 
we obtain 



W p {a) = 2 Z p (K,ix)J (2xs/a)xdx ) (7) 
Jo 
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It is important that the last equation has the partition 
function at imaginary value of the coupling constant. 
This should be understood as analytic continuation of 
Z p (K,x). 

There are several ways how one can compute W p (a). 
The most straightforward approach is to evaluate the co- 
efficients and thus determine all the moments of the 

(rt) 

distribution. Explicit expressions for Z^ can be ob- 
tained using orthogonal Jack-polynomials^. However, 
each of these coefficients is given as a series of products 
of T-functions and their evaluation becomes extremely 
cumbersome for n > 3. Another approach to finding 
W p (a) is to compute Z p (K,x) using the thermodynamic 
Bethe Ansatz for the quantum impurity problem 2 ^. This 
method works only for half integer K and requires so- 
lution of coupled integral equations. Besides, relating 
Z p (K, x) to the distribution function of interference am- 
plitudes requires analytic continuation of Z P (K, x) into 
the complex plane, x — > ix (see Eq. (0)), which introduces 
additional complications. In this paper we will use a dif- 
ferent method, which is based on studies of the integrable 
structure of conformal field theories24. In particular, it 
was shown that the vacuum expectation value of Bax- 
ter's Q operator, central to the integrable structure of 
the models, coincides with the grand partition function 
of interest up to an overall prefactor— (see supplemen- 
tary material for more details): 

Q™(c,X) = \^ K Z p (K,-ix), (8) 

where x is related to the spectral parameter A as x = 
7rA/ sin(7r/2i4T) and the central charge c = 1 — 6(2K — 
1/(2K)) 2 . It was conjectured in Refs. |25l26j that the 
vacuum expectation value Q p ac (X) is proportional to the 
spectral determinant of the single particle Schrodinger 
equation 

-fi£*(a:) + (x AK - 2 + ^4r^) = E*(x), (9) 

where I = ipK - \. So, Q vac (A) = X 4p / K D(pX 2 ), where 
p = (AK) 2 -^ K T 2 (1 - l/2iQ, D ( E ) is the spectral de- 
terminant denned as D(E) = Il^Li(l — E/E n ), and E n 
are the eigenvalues of 0. Thus, we get 

Z p (K,ix) = f[(l-^f\. (10) 

n— 1 ^ n ' 

To evaluate the distribution function we solve the 
Schrodinger Eq.® numerically. Details of the analyti- 
cal treatment and comparison with numerics will be re- 
ported elsewhere 2 ^. We checked the accuracy of the nu- 
merics as well as conjecture IjlOU by comparing coeffi- 
cients zf(K) and zf(K) evaluated for various K us- 
ing (i) the spectral determinant (e.g. Z% ~ ^2l/E n , 

z 4 } ~ dZ 1 / E n) 2 - E 1 /^) ). and (ii) the exact ex- 
pressions of Ref. [23 based on Jack polynomials. Wc 
found perfect agreement between the two methods. 



To compare distributions at different K with each 
other it is convenient to use a normalized interference 
amplitude a = a/ (a) = A 2 /(A 2 ) instead a. This change 
of variable is also convenient for comparison with ex- 
periments. The distribution function Wo (a) is shown in 
Fig- for several values of K. For K close to 1 (Tonks- 
Girardeau limit) Wo is a wide Poissonian function, which 
gradually narrows as K increases, finally becoming a nar- 
row 5-function at K — > 00 (the limit of noninteracting 
bosons). Interestingly, the distribution function remains 
asymmetric for arbitrarily large K. In fact we find that 
Wo(a — 1) tends to a universal scaling form, parameter- 
ized by a single number characterizing the width of the 
distribution: 5a = ■J (a 2 ) — 1 w n/\f&K (see Ref. 0). 
We conjecture that this limiting form of Wo is the Gum- 
bel distribution 28 , which frequently appears in problems 
of extreme value statistics^: Wo (a — 1) — » Wcia, K, 7), 
where 7 ~ 0.577 is the Euler gamma-constant and 

Wg(x, a, b) = aexp (ax — b — exp(as — b)) . (11) 

We plot the scaled distribution functions: 
SaW ((a — 1) /Sa). Note that Wo was multiplied 
by Sa to preserve the normalization condition (so that 
the total probability is equal to unity) for K = 10, 20, 40. 
For comparison in Fig. [3] we also present the scaled 
Gumbcl distribution. One can see that as K increases 
the function Wo indeed approaches to Wg- Gumbcl 
distributions are frequently associated with random 
walks in strongly correlated systems^. Indeed one 
can view the interference signal in Eq. P as a sum of 
contributions coming from different points along the 
condensates. For weak interactions (large K) these 
contributions are strongly correlated because the phases 
of each of the condensates only weakly fluctuate along 
z. Thus there is no surprise that Wo approaches the 
Gumbcl distribution. One can also understand this 
result noting that for K ^> 1 the distribution function of 
the interference amplitude is dominated by rare events 
which reduce the contrast. The Gumbel distribution 
was introduced precisely to describe rare events such as 
stock market crashes or earthquakes. In supplementary 
materials (sec. IV B|) we also discuss distribution func- 
tions for finite values of the observation angle (i.e. finite 

Interestingly, the distribution function W p (a) provides 
a very simple and convenient framework for describing 
both the partition function (0} and the expectation val- 
ues of the Q-operator. Indeed, W p (a) is a smooth well 
behaved function at all values of K and it can be easily 
approximated by simple analytic expressions. 

III. QUANTUM SIMULATION 

As we mentioned earlier, the distribution function 
W p (a) can be used to obtain the partition function 
Z p (K, x) (see Eq. @) describing a range of various prob- 
lems like quantum impurity in a one-dimensional elec- 
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FIG. 2: Evolution of the distribution function Wo (a) for dif- 
ferent values of K at p = 0. At larger values of K the function 
W p tends to the delta-function (see text for details) . 
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FIG. 3: Scaled distribution function SaWo ((& — 1)/Sa), 
where 5a is the width of the distribution, for large K. The 
function Wo is multiplied by 5a to preserve the total proba- 
bility, which must be equal to unity. The dashed and dotted 
lines correspond to different values of K. The solid line cor- 
responds to the conjectured Gumbel distribution. 



tron liquid, asymmetric Kondo problem, and dissipative 
tunneling. It is easy to see that momentum p in the 
interference experiments corresponds to the external ap- 
plied voltage in the impurity problem (ip ~ V). This 
follows from the interchanged roles of space and time in 
the two problems. Thus measuring W(a) experimentally 
and taking its integral transform one directly simulates 
these problems in or out of equilibrium. Moreover, af- 
ter substituting Io{2xyfa) — > Jo(2x-\/a) in Eq. @ one 



obtains the partition functions of the above models with 
imaginary coupling constants. Such models have been ac- 
tively investigated recently in the context of theories with 
PT-symmetric rather than Hermitean Hamiltonians^S. 

The partition function with the imaginary coupling 
also gives the expectation values of the Baxter Q- 
operator (see Eq. JSJl) corresponding to various con- 
formal field theories (CFTs). In general, such theo- 
ries are particularly important because many models of 
two-dimensional statistical mechanics, field theory, and 
many-body quantum systems at critical points can be de- 
scribed by some continuum theories having a property of 
conformal invariance. This leads to a description of crit- 
ical systems on the basis of conformal field theory, which 
basic ingredients are the central charge and conformal di- 
mensions. This set of data classify different universality 
classes, which often describe very different physical mod- 
els. The property of positivity of the central charge leads 
to a unitary theory. Physically, the central charge deter- 
mines the vacuum (Casimir) energy of the system and 
governs the finite-size scaling effects. On the other hand 
there is a class of models whose universality classes of 
critical behavior are described by the conformal-invariant 
models with negative central charges. These theories are 
nonunitary and their properties are thus very different 
from those described by positive central charges. 

Conformal field theories corresponding to Eq. JSJ are 
characterized by negative central charge c = 1 — 6(2K — 
1/(2K)) 2 and the highest weight A = (2piT) 2 + (c-l)/24 
(see Ref. |24|). These relations give c < -2 for K > 1. 
Theories with negative central charges appear in differ- 
ent contexts of statistical mechanics, stochastic growth 
models, 2D quantum gravity, models of 2D turbulence 
and even high-energy QCD. In particular, c = —2 CFT 
has the field-theoretical representation in terms of the 
ghost (anticommuting) fields and also corresponds to the 
critical behavior of the non-intersecting loop model on a 
2D lattice^ as well as to the special case of the stochas- 
tic Loewner evolution equation, describing the growth of 
random fractal stochastic conformal-invariant interface 
(see e.g. Ref. |23)- The classic example of CFT with neg- 
ative c is the Yang-Lee singularity c — —22/5 describing 
the critical behavior of the Ising model in imaginary mag- 
netic field. Models of 2D quantum gravity described by 
the fluctuating lattice geometry are related to negative 
c as well. Possibly, the high-energy limit of multicolor 
QCD is described by the integrable CFT with negative 
(or zero) central charge^. We illustrate the dependence 
of c on K as well as particular examples of models corre- 
sponding to different values of c in Fig. 21 The spectrum 
of Q-operator can be used to reconstruct the transfer 
matrices of the above mentioned negative c models. Par- 
ticular example of such procedure for c = — 2 universal- 
ity class was explicitly constructed in Ref. (24[. In this 
case only the vacuum expectation value was needed to 
reconstruct the whole transfer matrix. The transfer ma- 
trices contain all the information about the properties 
of underlying models. In this sense interference experi- 
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FIG. 4: Dependence of the central charge of conformal field 
theories corresponding to interference experiments on the 
Luttinger liquid parameter K. The red circles correspond 
to some particular examples discussed in the text. 

merits simulate these models. Experimentally, the central 
charge can be tuned by varying the interaction K and the 
scaling dimensions A of corresponding physical operators 
can be manipulated by changing the observation angle p. 
An interesting challenge here is the experimental deter- 
mination of nonvacuum values of the Q and T operators. 
This is an open question. 

Needless to say that the range of models mentioned 
above, which belong to nonunitary universality classes 
is difficult (if possible at all) to realize by other ways. 
The interference of condensates provides a possible and a 
plausible way to explore the interesting physics of various 
models ranging from statistical to high energy physics. 



IV. SUMMARY 

To summarize, in this paper we analyzed interference 
experiments between two independent one dimensional 
quasi-condensates. We computed the distribution func- 
tion of the amplitude of interference fringes relating this 
problem to the properties of Q operators of conformal 
field theories with negative central charges. We showed 
how one can use the distribution function of the inter- 
ference amplitude to reconstruct the partition function 
of a two-dimensional Coulomb gas confined to a circle. 
This partition function is related to a variety of statis- 
tical and field-theoretical models. Thus studying the in- 
terference distribution function experimentally, one can 
directly simulate these interesting models. 

We considered only a particular example of inter- 
ference between two one-dimensional condensates and 
showed the connection between the distribution function 
of fringe amplitudes and the properties of various mod- 
els. This analysis can be extended to other systems with 
quasi long-range order, e.g. to two-dimensional Bose sys- 
tems at finite temperature. One can expect that there 
will be analogous connections to different classes of prob- 



lems, some of which might not be exactly solvable. The 
interference experiments open new ways of solving these 
problems by direct simulation of the underlying models. 
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V. SUPPLEMENTARY MATERIAL 

A. Mathematical details 

Function Z(K,g) can be computed using the Thermo- 
dynamic Bethe Ansatz. This was used in Refs. [T5II22] to 
evaluate the current through a boundary impurity. How- 
ever to compute the distribution function W we need 
analytic continuation of the partition function Z(K,ig). 
Expression for Z(K, g) is given as a solution of coupled 
integral equations and performing its analytic continu- 
ation is not easy. In this paper we use an alternative 
approach. 

In a recent series of papers Bazhanov, Lukyanov, and 
Zamolodchikov explored an integrable structure of con- 
formal field theories focusing on connections to solvable 
problems on lattices^. Key ingredients of solvability of 
lattice models are is the so-called transfer matrix T(A) 
operators. These operators contain information about 
all integrals of motion as well as excitation spectra of 
the system. Transfer matrices are defined as a function 
of the so-called spectral parameter A (in the continuum 
limit A corresponds to rapidity) and commute for differ- 
ent values of A. The latter property is a direct manifes- 
tation of the existence of infinite number of commuting 
integrals of motion. In his studies of 8-vertex model, 
Baxter"^ introduced the operator Q(A) which helps to 
find an eigenvalues of T. Operators T and Q± satisfy a 
set of commutation relations, in particulars* 

T(A)Q±(A) =Q±(<zA) + Q±( ( T 1 A) (12) 

where q = exp(i7r/2if ). So T matrices can be obtained 
explicitly when one knows the Q operators. 

Operators A±(A) = Q±(A)A T4p / x act in the repre- 
sentation space of Virasoro algebra, which can be con- 
structed from the Fock space of bosonic operators a± n 
satisfying a n \p) = 0, (n > 0). The Fock vacuum state \p) 
is an eigenstate of the momentum operator, P\p) = p\p). 
For p = N/2 (N = 0,1,2...) the vacuum eigenvalues of 

the operator Aj-(A), A±(\)\p) = A^^lp), are given by 
( below we consider only the quantities with the + sub- 
script which correspond to the positive p) 

A {vac \\) = Z p (K,-ix) (13) 
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where 



A = — sin I I 



(14) 



7T \2K J 

The function A^ vac ^ has known large-A asymptotical 

\og{A vac ){\) ~ M{K) (-A 2 ) , (15) 
where the constant M(K) is given by 

v^r(^)[r(2§^!)]^ 

M(K) = UK - 2 ' 1 K 2K JL . (i 6 ) 

cos( 37 ^)r(^ T ) 

The function A^^iX) is entire function for K > 1 and 
is completely determined by its zeros A&, k — 0,1,.... 
Therefore A^ vac \\) can be represented by the convergent 
product 



A {vac ){\) = ]~[ (l- 



fc=0 



A 2 fe 



A {vac \0) = l (17) 



On the basis of analysis of a certain class of exactly 
solvable model, corresponding to the integrable pertur- 
bation of the conformal field theory, it was conjectured 
in24 that the so-called F-system and related T system 
(where Y = e <Lr and e r are the Bethe-ansatz energies 
parametrized by r, the nodes of the Dynkin diagrams) 
satisfy the same functional equations and possess the 
same analytical structure and asymptotics as the spectral 
determinant of the one-dimensional anharmonic oscilla- 
tor. Further, the same functional equations, analytical 
properties l|17f) and asymptotics l|15H are satisfied for the 
vacuum eigenvalues of Q-operator for special values of p 
and the latter are given by the spectral determinant of 
the following Schrodinger equation 



{~d 2 x + x 2a )^>{x) =EV(x). 
The spectral determinant is defined as 



^)=n(i-# 



E,, 



(18) 



(19) 



Soon after, in Ref. [26j, this conjecture has been extended 
to all values of p : 



A^ vac \\,p) = D( P X 2 ), 



(20) 



where now D(E) is the spectral determinant of Eq. @ 
with I = 4pK - 1/2. Here p = (4K) 2 - 1 ' ' K '[I '(1 - 
l/{2K))f. 

Typically, for n > 5 — 10 the spectrum of the equa- 
tion © is very well approximated by the standard WKB 
expression 



E n = e(K)(n-T l (K))-K- 



(21) 



where n = 1,2, ... Here ji(K) is the Maslov index. For 
I _ I 

4 2 ' 



1/2 < K < oo, ji(K) = i-L for K = oo, lt {K) = -1/2 



and for < K < \, 7/ (iT) = Note that the 

Eq. JSJ) has an interesting duality symmetry (generalizing 
the Coulomb-harmonic oscillator duality) which allows to 
relate the K > 1/2 and < K < 1/2 sectors. The point 
K = 1/2 is a self-dual point of this transformation. The 
function e(K) in Eq. I|21|l reads 



e(K) 



20FF(§ + jj^] 



r(i 



4K-2 



(22) 



In principle, the function ji(K) can be a smooth func- 
tion interpolating between limiting values given above 
and can be considered as a noninteger Maslov inde:»2&. 
This interpolation allows to (approximately) evaluate the 
partition function and the distribution function in many 
cases. In the limiting cases K — » 1 and K — >■ oo the WKB 
approximation gives the exact spectrum, which was dis- 
cussed in the main text. We point however, that the 
using the approximate WKB spectrum can result in non- 
physical results, e.g. negative values of the distribution 
function W(a), and thus have to be used with care. 



B. Analysis of distribution functions 

Analytic expressions for Z can be obtained explicitly 
for K — 1 since in this case equation © corresponds to 
a singular harmonic oscillator (with a parabolic poten- 
tial replaced by the potential x 2 + 1(1 + l)/x 2 , singular 
at x — 0). Eigenvalues of this Schrodinger equation are 
E n = An + 21 — 1, n = 1, 2, .... Exact result for the spec- 
tral determinant can be computed using the Weierstrass 
representation of the gamma function. We find, 

Z p (l,ix) = T(i + 2p)e~^ 2 [r(i + 2p - (23) 

where C is a nonuniversal constant which arises due to 
the logarithmic divergence of the integrals in Eq. © and 
involves the cutoff dependence of (a) at K = 1 . If we are 
interested in distances much larger than the cutoff then 
C> 1 and the function Z p becomes a simple gaussian. 
After the integral transform this leads to the Poissonian 
distribution W p (a). The case K — > oo is recovered by 
the infinite well potential for which the cigcncncrgies are 
given by zeroes of the Bessel function. Therefore 



Z p (K,ix) = T{ApK + l) X - 4pK .h pK {2a 



(24) 



From equations and l|24|) one sees that when p = 0, 
the distribution function is a delta function. As p in- 
creases, W p (a) rapidly broadens. In the limit of large 
K, the function W p depends only on the product Kp 
and takes a simple form: W p (a) « AKp (1 — a) AKp ^ 1 for 
a < 1 and zero for a > 1. This function is peaked at 
a = 1 for Kp < 1/4, it becomes a step function exactly 
at Kp — 1/4, and for Kp > 1/4 W p is a monotoni- 
cally decreasing function of a. When the product Kp 
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becomes large Kp > 1 the function W p becomes Poisso- 
nian: W p (a) w 4_Kpexp(— AKPa). In general, the ten- 
dency of broadening of the distribution function remains 
true for all values of K . In particular, we show the be- 
havior of W at fixed K = 2 for various p in Fig. JSJ. 



FIG. 5: Evolution of distribution function at K = 2 for dif- 
ferent values of index p (i.e. for different observation angles). 
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